DISCRETE PERIOD MATRICES AND RELATED 

TOPICS 



CHRISTIAN MERCAT 

Abstract. We continue our investigation of Discrete Riemann 
Surfaces with the discussion of the discrete analogs of period ma- 
trices, Riemann's bihnear relations, exponential of constant argu- 
ment, series and electrical moves. We show that given a refining 
sequence of critical maps, the discrete period matrix converges to 
the continuous one. 



1. Introduction 

The notion of discrete Riemann surfaces was defined in The 
interesting paper |^ initiated a renewed interest in the subject. In 
their paper, R. Costa-Santos and B. McCoy observed numerically that 
certain pfaffians intervening in a dimer or critical Ising model converge 
at the thermodynamic limit to a certain (power of a) theta function at 
the origin. They computed the period matrix needed to define the theta 
function using discrete holomorphy method. The present paper aims 
at putting their work in a more general theoretical framework. Most 
of the results in this paper are a straightforward application of the 
continuous theory |^ together with the results in |]I|, 0, |[, to which 
we refer for details. We define the discrete period matrix, which is 
twice as large as in the continuous case: the periods of a holomorphic 
form on the graph and on its dual are in general different, but the 
continuous limit theorem, given a refining sequence of critical maps, 
ensures that they converge to the same value. The main tool is the 
same as in the continuous case, the Riemann bilinear relations. We 
define the discrete exponential of a constant argument on a critical 
map and explore its properties as well as its link with series. The tools 
for that purpose, which needs more investigations in its own right, are 
the electrical moves. 



2. Discrete Riemann surfaces 

We recall in this section basic definitions and results from where 
the notion of discrete Riemann surfaces was defined. We are interested 
in discrete surfaces given by a cellular decomposition <) of dimension 
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two, where all faces are quadrilaterals (a quad-graph [0]). It defines, 
up to homotopy and away from the boundary, two dual cellular de- 
compositions r and r*. Edges Tl are dual to edges Fi, faces T2 are 




0. The vertex dual to a face. 1. Dual edges. 2. The face dual to a vertex. 



Figure 1. Duality. 

dual to vertices Fq and vice-versa. Their union is denoted the double 
A = F U F*. A discrete conformal structure on A is a real positive 
function p on the unoriented edges satisfying p(e*) = l/p(e). It defines 
a genuine Riemann surface structure on the discrete surface: Choose 
a length 6 and realize each quadrilateral by a lozenge whose diagonals 
have a length ratio given by p. Gluing them together provides a fiat 
riemannian metric with conic singularities at the vertices, hence a con- 
formal structure [§]. It leads to a straightforward discrete version of 
the Cauchy- Riemann equation. A function on the vertices is discrete 
holomorphic iff for every quadrilateral {x , y , x' , y') G <C>2, 

(2.1) f{y') - f{y) = tp{x,x') {f{x') - f{x)) . 



We recall elements of de-Rham cohomology, doubled in our context: 
The complex of chains C(A) = Co (A) © Ci(A) © C2(A) is the vec- 
tor space span by vertices, edges and faces. It is equipped with a 
boundary operator d : Cfc(A) Cfc-i(A), null on vertices and fulfilling 
= 0. The kernel ker d =: Z,{A) of the boundary operator are the 
closed chains or cycles. Its image are the exact chains. It provides the 
dual spaces of forms, called cochains, C^{A) := Hom(Cfc(A),C) with a 
coboundary d : C'^(A) — > C^^^(A) defined by Stokes formula: 

J df:=f {d{x, x')) = fix') - fix), jj da:= ja. 

{x,x') F dF 

A cocycle is a closed cochain and we note a G Z'^iA). 
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Figure 2. The discrete Cauchy-Riemann equation. 



Duality of complexes allows us to define a Hodge operator * on forms 



by 



(2.2) 



* : C\A) 
C°(A) 3 f 

C\A) 3 a 

C\A) 3 00 



*/ := fiF* 
-p{e 
*UJ : {*uj){x) := 



F 



a, 



It fulfills *^ = (— Id)^ and defines A := —d *d* — *d*d, the usual 
discrete Laplacian. Its kernel are the harmonic forms. The discrete 
holomorphic forms are special harmonic forms: a 1-form 



(2.3) 



a G C^(A) is holomorphic iff rfa = and * a = —ia, 



that is to say if it is closed and of type (1, 0). We will note a G i7^(A). 
A function / : Aq — C is holomorphic iff df is holomorphic and we 
note / e Q\A). 

In the compact case, the Hodge theorem orthogonally decomposes 
forms into exact, coexact and harmonic; harmonic forms are the closed 
and co-closed ones; and harmonic 1-form are the orthogonal sum of 
holomorphic and anti-holomorphic ones. We studied discrete Hodge 
theory for higher dimensional complexes in 

We construct a wedge product on (} such that c?^ is a derivation 
for this product A : C"=«» x C\0) It is defined by the 
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following formulae, for f,ge C^O), a,(3e C\<» and uj G C^O): 
if ■ 9)ix) ■.=f{x) ■ g{x) for X G Oo, 

/ ■ a ■= ^^^^ ^ ^^y^ I {x,y)a for (a;, G Oi, 

{x,y) ^ J 

aAp:=^^ j a j (3- j a j (3 

(X1,X2,X3,X4) ^~^iXk-l,Xk) (Xk,Xk + l) {Xk + l,Xk) {Xk,Xk-l) 

f{xi)+f{x2)+f{x3)+fix4) f f 

{xi,X2,X3,X4,) (xi,X2,X3,X4) 

for (Xi,X2,X3,X4) G <^2- 

A form on <) can be averaged into a form on A: This map A from 
C*(<^) to C*(A) is the identity for functions and defined by the following 
formulae for 1 and 2-forms: 

. . . . . \ 

(2.4) / AK): ^ 




(2.5) ffM^0)--=li2 



2\J ' J ' J ' J I 

{x,x') \x,y) {y,x') {x,y') {y',x') j 

2^ // 

fc=i 

^* (xk,yk,x,yk-i) 

where notations are made clear in Fig. |^. The map A is neither injective 
nor surjective in the non simply-connected case, so we can neither 
define a Hodge star on <^ nor a wedge product on A. Its kernel is 
Ker (A) = Vect (c?oe), where e is the biconstant, +1 on F and —1 on 
F*. But d\A = Ad(^ so it carries cocycles on (} to cocycles on A and 
its image are these cocycles of A verifying that their holonomies along 
cycles of A only depend on their homology on the combinatorial surface: 
Given a 1-cocycle fi G Z^{A) with such a property, a corresponding 
1-cocycle u G Z^{(}) is built in the following way, choose an edge 
{xo,yo) G <^i; for an edge {x,y) G <^i with x and xq on the same leaf 
of A, choose two paths Xx,xo and Xy^^y on the double graph A, from x 
to Xq and yo to y respectively, and define 



(2.6) / z/:=/ fi- f 

where [7] = [A^.^^.^ + (xq, yo) + Xy^^y + (y, x)] is the class of the full cycle 
in the homology of the surface. Changing the base points change /i by 
a multiple of (i^e. 
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It follows in the compact case that the dimensions of the harmonic 
forms on (} (the kernel of AA) modulo de, as well as the harmonic forms 
on A with same holonomies on the graph and on its dual, are twice the 
genus of the surface, as expected. Unfortunately, the space Im A = 
?i-'-©Im d is not stable by the Hodge star *. We can nevertheless define 
holomorphic 1-forms on (} but their dimension can be much smaller 
than in the continuous, namely the genus of the surface. Criticality 
provides conditions which ensure that the space *Im A is "close" to 
Im A. 



X 



y 
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x' 



'■o' 

y 




Figure 3. Notations. 



We construct an heterogeneous wedge product for 1-forms: with 
a,pe C\A), define a A (3 e C\0) by 



(2.7) 



aAp:= J a j (3 + J a j (3. 

{x,y,x',y') {x,x') {y,y') {y,y') {x' ,x) 



It verifies A{a^) A A{i3^) = A (3^, the first wedge product being 
between 1-forms on A and the second between forms on (). The usual 
scalar product on compactly supported forms on A reads as expected: 



2.8) (a, (3) := ^ p(e) (^J^a^ (^j^(3^ = JJc^^*P 



3. Period matrix 



We use the convention of Farkas and Kra |Q , chapter III, to which we 
refer for details. Consider (<),p) a discrete compact Riemann surface. 
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3.1. Intersection number, on A and on (}. In [|T], for a given simple 
(real) cycle C G Zi{A), we constructed a harmonic 1-form rjc such that 
§^ rjc counts the algebraic number of times A contains an edge dual to 
an edge of C. It is the solution of a Neumann problem on the surface 
cut open along C. It follows from standard homology technique that 
r]c depends only on the homology class of C (all the cycles which differ 
from C by an exact cycle dA) and can be extended linearly to all cycles; 
it fulfills, for a closed form 6, 

(3.1) = // VcA0, 

Jc J J 

and a basis of the homology provides a dual basis of harmonic forms 
on A. Beware that if the cycle C G ZiiT) is purely on F, then this 
form r^dp = is null on F. 

The intersection number between two cycles A, B E ^i(A) is defined 

as 

(3.2) A-B:= JJ^7]aAVb. 

It is obviously linear and antisymmetric, it is an integer number for 
integer cycles. Let's stress again that the intersection of a cycle on F 
with another cycle on F is always null. A cycle C G Zi{(}) defines 
a pair of cycles on each graph Cr G Zi(T), Cr* G Zi(T*) which are 
homologous to C on the surface, composed of portions of the boundary 
of the faces on A dual to the vertices of C. They are uniquely defined 
if we require that they lie "to the left" of C as shown in Fig.^. By 
the procedure ( ^.6] ) applied to rjc^ + rjc^^ , we construct a 1-cocycle 
rjc G -Z'^(<)) unique up to de, and since V^, de f\6 = Eq. ( |3.2| ) defines 
an intersection number on Zi{()). Unlike the intersection number on A, 
this one has all the usual expected properties. In particular Eq. ( |3.2D 
holds for G Zi(^). 

3.2. Canonical dissection, fundamental polygon. The complex 
<C> being connected, consider a maximal tree T C <0i, that is to say T 
is a Z2-homologically trivial chain and every edge added to T forms a 
cycle. A canonical dissection or cut-system K of the genus g discrete 
Riemann surface <) is given by a set of oriented edges {ek)i<k<2g such 
that the cycles K C (TUe/c) form a basis of the homology group Hi{(}) 
verifying, for 1 < k,i < g 



(3.3) Kfe-K, = 0, Kfc+3-K,+3 = 0, ■ ^e+g = Sk,e. 
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Figure 4. A path C on () defines a pair of patlis Cr 
and Cr* on its left. 

Tliey actually form a basis of the fundamental group 7ri(<^) and the 
defining relation among them is (noted multiplicatively) 

9 

(3.4) n^'=^'=+A"'^^"+.=i- 

k=l 

The construction of such a basis is standard and we won't repeat the 
procedure. What is less standard is the interpretation of Eq. ( |3.4D in 
terms of the boundary of a fundamental domain, discretization intro- 
duces some subtleties. 

Considering T U 6fc clS cl rooted graph, we can prune it of all its 
pending branches, leaving a simple closed loop K^, attached to the 
origin O by a simple path (see Fig. ^), yielding the cycle N^. These 
three cycles are deformation retract of one another, C C T U 
hence are equal in homology. 

In the continuous basis of the homology can be realized by 

2g simple arcs, transverse to one another and meeting only at the base 
point. It defines an isometric model of the surface as a fundamental 
domain homeomorphic to a disc and bordered by Ag arcs to identify 
pairwise. In the discrete case, by definition, the set of the cellular 
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Figure 5. 1. A maximal rooted tree in a quadrilateral 
decomposition of the torus. 2. An additional edge defines 
a rooted cycle Ki, pruned of its dangling trees. 3. Its un- 
rooted version, the simple loop . 



complex minus the edges taking part into the cycles basis is homeomor- 
phic to a disc hence the surface is realized as a polygonal fundamental 
domain M. whose boundary edges are identified pairwise. 

But it is sometimes impossible to choose a basis of the homology 
verifying ( p.3|) by simple discrete cycles which are transverse to one 
another. If the path is not empty, the cycle is not even simple. 
Moreover, some edges may belong to several cycles. In this case, the 
edges on the boundary of this fundamental polygon can not be assigned 
a unique element of the basis or its inverse, and therefore can not be 
grouped into only Ag continuous paths to identify pairwise but more 
than Ag. 

In fact, the information contained into the basis K is more than 
simply this polygon, the set of edges composing the concatenated cycle 

(3.5) (Ki, Vi, N2, . . . , K^-i) 

encodes a cellular complex M.^ which is not a combinatorial surface 
and consists of the fundamental polygon M. plus some dangling trees, 
corresponding to the edges which belong to more than one cycle or 
participate more than once in a cycle (the paths Afc), as exemplified 
in Fig.p. By construction, the edge Ck belongs to the cycle only, 
hence these trees are in fact without branches, simple paths whose 
only leaf is the base point O. To retrieve the surface, the edges of this 
structure are identified group-wise, an edge participating k times 
in cycles will have [fc/2] -|-2 representatives to identify together, two on 
the fundamental polygon and the rest as edges of dangling trees. 

Ehminating repetition, that is to say looking at ( p3D not 
quence of edges but as a simplified cycle (or a simplified word in edges), 
thins Ai+ into A4, pruning away the dangling paths. The fundamental 
polygon boundary loses its structure as Ag arcs to be identified pairwise, 
in general a basis cycle will be disconnected around the fundamental 
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Figure 6. Three different fundamental polygons of a 
decomposition of the torus {g = 1) hj three quadrilat- 
erals: 1. The standard fundamental domain where the 
ig paths are not adapted to <C>. 2. A^_|_ is composed of 
edges of <^ composing 4g arcs (which may have portions 
in common) to identify pairwise, each edge corresponds 
to an element of the basis K or its inverse, except for 
edges of "dangling trees" which are associated with two 
such elements. 3. A4 is composed of edges of <) compos- 
ing more than ig arcs to identify pairwise, there is no 
correspondence with a basis of cycles. 



domain and a given edge can not be assigned to a particular cycle. 
This peculiarity gives a more complex yet well defined meaning to the 
contour integral formula for a 1-form 9 defined on the boundary edges 
of M+, 

(3.6) / e = y^l e + 

JdM J-^k 

This basis gives rise to cycles and '^^* whose homology classes form 
a basis of the group for each respective graph, that we compose into 
defined by 

(3.7) = C <, = C 

for 1 < k < g so that while the intersection numbers matrix on <0 is 
given by the 2g x 2g matrix 

(3.8) (^.•^.W=(_°j o)' 
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the intersection numbers matrix on A is the Ag x Ag matrix with the 
same structure 



(3.9) 



r 




r* 


r 




/o 





/ 


o\ 


r 











/ 


r* 


-/ 











r* 


^0 


-/ 





V 


r 



3.3. Basis of harmonic forms, basis of holomorphic forms. We 

define a^, the basis of real harmonic 1-forms, dual to the homology 
basis K^, 



at 



(3.10) 

which verify 



(3.11) 



a 



k+2g 



and 

for 1 < A; < 2^ 



A 
k+2g 



'kl, 



and dually, the intersection matrix elements are given by 



(3.12) 



^k ■ 



* Of 



On (}, the elements 
defined up to de, verify A{al 





Vn,+^ and a^^g 



for 1 < /c < 5f, 



a 



and form a basis of the cohomology on dual to K as well, 



k+2g 



+ a 



k+3g 



(3.13) 



a 



k+g 



and 

for 1 < k < g, 



they fulfill the first identity in Eq.( p.l2| ) but the second is meaningless 
in general since * can not be defined on <^. We will drop the mention 
A when no confusion is possible. 



Proposition 3.1. The matrix of inner products on A, 
(3.14) 

ak A*ai= < 7^*+2« 
is a real symmetric positive definite matrix. 



[ak, ai)k,£ 



A D 
B C 
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Proof It is real because the forms are real, and symmet- 

ric because the scalar product (|2.8| ) is skew symmetric. Definition 
Eq. ( |3.10|) and Eq. lead to the integral formulae. Positivity fol- 
lows from the bilinear relation Eq. ( [4.2| ): for 6 = Ylk=i^k ctki with 

a e C, E?=i lap > 0, 




(3.15) = J2^k^iiak,ai) >0. 





The form ak is supported by only one of the two graphs F or T*, 
the form *ak is supported by the other one, and the wedge product 
6'r A 6*1. = is null for two 1-forms supported by the same graph. 
Therefore the matrices A and C are g x gi-block diagonal and B is 
ant i- diagonal. 



(3.16) 
A = 



At 







Ar^ 



B 





Btt* 







C 











The matrices of intersection numbers (|3.9| ) and of inner products differ 
only by the Hodge star *. Because * preserves harmonic forms and the 
inner product, we get its matrix representation in the basis a. 



(3.17) 



-D A 
-C B 



and because *^ = — 1, 



(3.18) 
(3.19) 
(3.20) 



B^ -C-A + I 
A-B 
C-'B 





'B-A 
B-C. 



On <0, while the Hodge star * can not be defined, we can obviously 
consider the following positive scalar product on the classes of closed 
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forms modulo de, to which the set (a^) belong: 



[a 



{A{aO),A{f30)) 



t 



E 



p=p{x,x'),p*=p(y,y') 



(a; ,y ) 



a 



a 



I{x',y') ^ 



^ +P+P* 


+P-P* 


-p-p* 


-p+p*^ 




+p+p' 


-p+p* 


* 

-p-p 


* 

-p-p 


-p+p* 


+P+P* 


+P-P 


y-p+p* 


-p-p* 


+P-P* 


+p+p*y 



and it yields 

(3.21) («^«f),,, 



which, in general, can not be understood as the periods of a set of 
forms on <0> along the basis K. 

Let's decompose the space of harmonic forms into two orthogonal 
supplements, 

(3.22) n\K)='H\®^n\ 

where the first vector space are the harmonic forms whose holonomies 
on one graph are equal to their holonomies on the dual, that is to say 

(3.23) n\ := Vect (a^ + ak+g, I < k < g or 2g < k < 3g). 
Definition ( p. 101 ) and Eq. ( p.lD imply that 

(3.24) n\ = Vect {*ak - *ak+g, l<k<goi2g<k< 3g). 

These elements in the basis (a^ + a^+g, ; ctk — c^k+g) ioi 1 < k < g and 
2g < k < 3g, are represented by the following invertible matrix: 

^ I — ^ By*y — Ay* ^ 

/ Cp — Cp* Byy* — By*y 

*5rr- + ^Byt Ay + Ay* 

^0 Cy + Cy* Byy* ~\~ By*y J 

It implies in particular that the lower right g x g block is invertible, 
therefore so is Eq. ( p.21|) . 



(3.25) 



3.4. Period matrix. 

Proposition 3.2. The matrix U ■ 
of the basis of holomorphic forms 



C ^ ■ {i — B) is the period matrix 



29 



(3.26) 



5Z '^k} 



kl Oii+2g 



I{y ,x') 
I(x',y') 
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in the canonical dissection that is to say 

Jn^ [ilk-2g/ for 2g <k< Ag, 
and n is symmetric, with a positive definite imaginary part. 

The proof is essentially the same as in the continuous case |Q and 
we include it for completeness. 

Proof ^7^ . Let Uj := aj + i* aj for 1 < j < 4(7. These holomorphic 
forms fulfill 

(3.28) Pkj := ^{uJk, uJj) = (ttfc, ttj) + i ("fe, - * "i) 



(3.29) 



-I 



L ^Wfc, l<J<2g, 



' ^fc' 2^ < J < Ag. 



P is the period matrix of the forms [u) in the homology basis K. The 
first 2g forms {uJj)i<j<2g are a basis of holomorphic forms. It has the 
right dimension and they are linearly independent: 

2a 2g / 2g 2g \ 

^{Xj + ilj,j){aj + i * aj) = I (Aj + /^fc %fc) aj + A^fc Cj,fc a2g+i j 

j=l j=l \ k=l k=l / 

2g / 2g 2g \ 

(3.30) +i I + ^^k Bj^k) aj + ^Xk Cj^k a2g+j j 

j=l \ k=l k=l J 

is null, for A, /i G R only when A = /i = because C is positive 
definite. Similarly for the last 2g forms. The change of basis iC^^ on 
them provides the basis of holomorphic forms {(). The last 2g rows of 
P is the 2g x 4g matrix {B — i I, C) hence the periods of {() in K are 
given by (/,n). 

The first identity in Eq.( p.27|) uniquely defines the basis ( and a 
holomorphic 1-form is completely determined by whether its periods 
on the first 2g cycles of K, or their real parts on the whole set. 

Notice that because C is g x g block diagonal and B is ant i- diagonal, 
n is decomposed into four g x g blocks, the two diagonal matrices form 
iC~^ and are pure imaginary, the other two form —C^^ ■ B and are 
real. 



(3.31) n 
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Therefore the holomorphic forms (k are real on one graph and pure 
imaginary on its dual, 

(3.32) i<k<g aGC^(r)©zC^(r*) 

g<k<2g ^ CkeC^ir)®tCl,iT). 

We will call 

(3.33) Hr = + Hi, 

the period matrix on the graph T the sum of the real periods of (k, 
^ ^ k < g, on r, with the associated pure imaginary periods on the 
dual r*, and similarly for (k, g < k < 2g, the period matrix on T*. 

It is natural to ask how close Ilr and Ilr* are from one another, and 
whether their mean can be given an interpretation. Criticality 0] 
redefined in Sec. |^, answers partially the issue: 

Theorem 3.1. In the genus one critical case, the period matrices Ur 
and Ilr* are equal to the period matrix lis of the underlying surface S. 
For higher genus, given a refining sequence Pk) of critical maps 
of S, the discrete period matrices Ilpfc and Ilp.fc converge to the period 
matrix lis • 



Proof \3l\ . The genus one case is postponed to Sec. p^ . The contin- 
uous limit comes from techniques in , to be developed in Q which 
prove that, given a refining sequence of critical maps, any holomorphic 
function can be approximated by a sequence of discrete holomorphic 
functions. Taking the real parts, this implies as well that any harmonic 
function can be approximated by discrete harmonic functions. In par- 
ticular, the discrete solutions fk to a Dirichlet or Neumann problem on 
a simply connected set converge to the continuous solution / because 
the latter can be approximated by discrete harmonic functions gk and 
the difference fk — gk being harmonic and small on the boundary, con- 
verge to zero. In particular, each form in the basis {af), provides a 
solution to the Neumann problem Eq. ( |3.13 ) and a similar procedure. 



detailed afterwards, define a converging sequence of forms (f, yielding 
the result. <^ 
We can try to replicate the work done on A on the graph <^. A 
problem is that Ar + A-p and Cr + Cr need not be positive definite. 
Moreover, the Hodge star * doesn't preserve the space {A{a^)) of har- 
monic forms with equal holonomies on the graph and on its dual, so 
we can not define the analogue of a + i * a on (}. We first investigate 
what happens when we can partially define these analogues: 

Assume that for 2g < k < 3g, the holonomies of *ak on T are equal 
to the holonomies of *ak+g on T*, that is to say Cr = Cr* ='■ and 
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_Dpp* = £)p*p =: ^D^. It implies that the transposes fulfill i?rr* 
_ well. ... ...... ........ 



-Br*r =: |-Bo as well. We can then define G Z^{(}) such that 



A{f^k-q) ~ *'^k+g, uniquely up to de. The last g columns of 



the matrix of scalar product Eq. (p.21|) are related to their periods in 



the homology basis in a way similar to Eq. ( p.l4| ). By the same 
reasoning as before, the forms 

(3.34) = ^2 C^\,e - > l<k<9 



verify A((^) = ^^'''^^ and have periods on given by the identity 
for the first g cycles and the following g x g matrix, mean of the period 
matrices on the graph and on its dual: 



(3.35) = C^\i- B. 



The same reasoning applies when the periods of the forms *ak on 
the graph and on its dual are not equal but close to one another. In the 
context of refining sequences, we said that the basis (af), converges to 
the continuous basis of harmonic forms defined by the same Neumann 
problem Eq. (|3.13|) . Therefore 



(3.36) Cp-Cp, =o(l), Brr* - Br^r = o{l). 

A harmonic form u^+g = o(l) on T* can be added to *ak+g such that 
there exists G Z^{(}) with A{l3f_g) = *ak+g + jyk+g, yielding forms 
verifying A{(^) = ^{(k + Cfc+g) + o{l) and whose period matrix is 
+ o(l). Since the periods of converge to the same periods as its 
continuous limit, this period matrix converges to the period matrix lis 
of the surface. Which is the claim of Th. 

In the paper R. Costa-Santos and B. McCoy define a period 
matrix on a special cellular decomposition F of a surface by squares. 
They don't consider the dual graph T*. Their period matrix is equal to 
one of the two diagonal blocks of the double period matrix we construct 
in this case. They don't have to consider the off-diagonal blocks because 
the problem is so symmetric that their period matrix is pure imaginary. 

3.5. Genus one. Criticality solves partially the problem of having two 
different g x g period matrices instead of one since they converge to 
one another in a refining sequence. However, on a genus one critical 
torus, the situation is simpler: The overall curvature is null and a 
critical map is everywhere fiat. Therefore the cellular decomposition is 
the quotient of a periodic cellular decomposition of the plane by two 
independant periods. They can be normalized to (1, r). The continuous 
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period matrix is the 1 x 1-matrix r. A basis of the two dimensional 
holomorphic 1-forms is given by the real and imaginary parts of dZ on 
r and r* respectively, and the reverse. The discrete period matrix is 

the 2x2 matrix ^ and the period matrices on the graph 

and on its dual are both equal to the continuous one. 

For illustration purposes, the whole construction, of a basis of har- 
monic forms, then projected onto a basis of holomorphic forms, yielding 
the period matrix, can be checked explicitely on the critical maps of the 
genus 1 torus decomposed by square or triangular /hexagonal lattices: 

Consider the critical square (rectangular) lattice decomposition of a 
torus <^ = (Ze*^ + Ze~*^)/(2pe*^ + 2g e~*^), with horizontal parameter 
p = tan6 and vertical parameter its inverse. Its modulus is r = ^e^*^. 
The two dual graphs F and F* are isomorphic. An explicit harmonic 
form a\ is given by the constant l/2p on horizontal and downwards 
edges of the graph F and on all the other edges. Its holonomies are 
1 and on the p, resp. q cycles. Considering l/2q and the dual graph, 
we construct in the same fashion a2,a\*,a2 ■ The matrix of inner 
products is 



(3.37) iak,ai)k,e 



I 2 COS 2^\ 



sin 29 



2 cos2e 

cos2^ £ 
1 



Vos29 2 J 

using = 1/ sin 26* and ^^^^ = —1/ tan 26* so that the period 

matrix is 

(3.38) ^^q^^sm29 cos 2^ 

P 



cos 2^ i sin 2^ 



Therefore there exists a holomorphic form which has the same periods 
on the graph and on its dual, it is the average of the two half forms 
of Eq. ( |3.27| ) and its periods are (1, |e^*^) along the p, resp. q cycles, 
yielding the continuous modulus. This holomorphic form is simply the 
normalized fundamental form 

In the critical triangular /hexagonal lattice, we just point out to the 
necessary check by concentrating on a tile of the torus, composed of two 
triangles, pointing up and down respectively. We show that there exists 
an explicit holomorphic form which has the same shift on the graph 
and on its dual, along this tile. Let P-,p\ and p/ the three parameters 
around a given triangle. Criticality occurs when p^ p\ + p\p/ + p/ p_ = 
1. The form which is 1 on the rightwards and South- West edges and 
elsewhere is harmonic on the triangular lattice. Its pure imaginary 
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companion on the dual hexagonal lattice exhibits a shift by i p\ in the 
horizontal direction and i {p\ + p_) in the North-East direction along 
the tile. Dually, on the hexagonal lattice, the form which is p\ p_ along 
the North-East and downwards edges and 1 — p\p- along the South- 
East edges, is a harmonic form. Its shift in the horizontal direction is 
1, in the North-East direction 0, and its pure imaginary companion on 
the triangular lattice exhibits a shift by i p\ in the horizontal direction 
and i {p\ + P-) in the North-East direction along the tile as before. 
Hence their sum is a holomorphic form with equal holonomies on the 
triangular and hexagonal graphs and the period matrix it computes is 
the same as the continuous one. This simply amounts to pointing out 
that the fundamental form dz can be explicitely expressed in terms of 
the discrete conformal data. 



4. Bilinear relations 



Proposition 4.1. Given a canonical dissection for two closed forms 

e,e'ez\o), 



j=l -'^j+a ''^3+a 



(4.1) jj A0' = J2yj> e (f> e' - f o f o' 

for two closed forms 6,6' G Z^{A), 

(4.2) // 6 A6' = y" ( { 6 { 6' - I 6 I 6' 



Proof |7T7| . Each side is bilinear and depends only on the cohomol- 
ogy classes of the forms. Decompose the forms onto the cohomology 
basis (afc). On A, use Eq (|3.12|) for the LHS and the duality property 
Eq. ( p. lip for the RHS. On ^, use their counterparts. <(> 

Notice that for a harmonic form 6 G 7i^(A), the form *6 is closed as 
well, therefore its norm is given by 

(4.3) 6 el-L\k) ^ \\6\\'^=y^\ I 6 i *d - (f 6 l *e 

j = l \J\^, J^j + 2a J^j + 2a J^j 



5. Criticality and exponentials 



5.1. Integration at criticality. We call a local map Z : (} D U ^ 
C from a connected, simply connected subcomplex to the euclidean 
plane, critical iff it is locally injective, orientation preserving and the 
faces of <) are mapped to rhombi such that the ratio of the diagonal 
euclidean lengths is given by p. The common length 6 of these rhombi 
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is a characteristic of the map. We showed in [|T], ^ that a converging 
sequence of discrete holomorphic forms, on a refining (the lengths 6 
go to zero and faces don't collapse) sequence of conformal maps of the 
same Riemann surface, converge to a genuine holomorphic form. The 
converse is also true, every holomorphic form on a Riemann surface 
can be approximated by a converging sequence of discrete holomorphic 
forms given a refining sequence of critical maps. And such a refining 
sequence exists for every Riemann surface. 

Of course, a very natural atlas of maps of a discrete Riemann surface 
is given by the applications from every face, considered as a subcom- 
plex, to a well shaped rhombus in C. Such an application is only well 
defined if two edges of the same face are not identified. They then 
form an atlas of the combinatorial surface and provide it with a gen- 
uine structure of Riemann surface, namely a riemannian fiat metric 
with conic singularities 0. A discrete Riemann surface is called criti- 
cal or critically flat if its conic singularities are all multiple of 27r. 

We've seen that given a function / and a 1-form dZ, the residue the- 
orem implies that if they are both holomorphic then f dZ is a closed 
1-form. The key point about criticality is that, given a critical map Z 
and a discrete holomorphic function /, one can construct the holomor- 
phic 1-form / dZ 



Reciprocally, every holomorphic 1-form is locally of such a form, unique 
up to e. The change of map acts as expected (see Prop. (|7.3|)) so for 
a critical discrete Riemann surface, where the conic singularity angles 
are multiple of 27r, / dZ can be continued as a well defined 1-form on 
the universal covering of (} wherever / is defined. Be careful that for 
another holomorphic form 6 G f^(0), the 1-form fOE Z^{(}) is only 
closed but not holomorphic in general. 

5.2. Exponential. As an interesting example of discrete holomorphic 
function we present here the discrete exponential of a constant argu- 
ment. It is the first step in trying to define the exponential map needed 
to provide Abel's correspondence between divisors and meromorphic 
functions. 

Different attempts have been made to define a discrete exponential, 
see [D, |lOl and references therein. 

Definition 5.1. Let Z a critical map and O G Aq such that Z{0) = 0, 
the origin. For A G C, define the holomorphic function exponential 



(5.1) 
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Exp{:\:) G ^^°(A), denoted z t— > Exp{:\: z), by 
jExp{:X:0) = 1 



(5.2) 



^dExp{:\: z) = \Exp{:\: z)dZ. 



Proposition 5.2. //|A| ^ 2/5, Exp{:\:) is a well defined holomorphic 
function. 



Proof ^7^ . If (x,?/) G <)i is an edge, Eq.(|5.2D reads 
(5.3) 

Exp(:A: ,) - Exp(:A: x) = A Exp(:A: + Exp(:A: .) ^^^^^ _ ^^^^^ 

SO that 

As the map Z is critical, on the face (x, y, x', y') G <C>2, Z{y) — Z{x) = 
Z{x') — Z{y') as well as Z{y') — Z{x) = Z{x') — Z{y) so that the product 
of the four terms along the face 

2 + \iy-x) 2 + \{x'-y) 2 + Xjy' - x') 2 + Xjx - y') ^ ^ 
^ ■ ^ 2~\{y-x) 2- \{x' -y)2- \{y' - x') 2- \{x - y') 

where we wrote Z{z) as z for readability. So Exp(:A:) is a well defined 
function on the simply connected sub-complex U C 0, its value at a 
point is the product of the contributions of each edge along a given 
path connecting it to the origin, and the result doesn't depend on the 
path. It is also holomorphic as 

Exp(:A: y ) - Exp(:A: y) 2-\{y'-x) 2-\(y-x) 



Exp(:A:a;') -Exp(:A:a;) ^+Ky-^) 2+Kx'-y) _ ^ 

(5.6) = — -=ip{x,x'). 



We note that it is logical to state e = Exp(:oo:) = ±1 on F and T* 
respectively. 

On the rectangular lattice (} = 5(e~*^Z + e*^Z), the exponential is 
given explicitly, for z = 5{ne~^^ + me*^) by 

(5.7) = exp(Az) + 0(5^) 
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as 6 goes to zero keeping z E C fixed, because (l+x)" = exp {n log(l + x)) = 
exp(nx) + O {{nxY). Similar results are obtained for other lattices, see 
Fig. 1^ for the comparison of the discrete and continuous exponentials 
on the triangular /hexagonal double. The discrepancy between the two 
increases with the modulus of the parameter and decreases as the mesh 
is refined. More generally, we will prove in ||^ the same behavior for 
any refining sequence of critical maps of a critical discrete Riemann 
surface. 

5.3. Change of base point. Given another critical map ( : V ^ C 
with an origin b G Aq, there exists a complex number a such that 
( = a{Z — 6) on a connected component of the intersection U (IV. If 
the discrete Riemann surface is itself critical and if ?7 U is contained 
in a ball then this number a is the same for all connected components. 
As d( = a dZ, it follows that 



as expected. Hence, on a critical discrete Riemann surface, the expo- 
nential can be continued as a well defined holomorphic function on the 
universal covering of the discrete surface. 

5.4. Train-tracks. For A with a norm equal to 2/5 in a given critical 
map, the associated exponential may not be defined as straightfor- 
wardly. We prove that in the compact case, in general it can only be 
defined as zero: 

We name the equivalence class of an (oriented) edge, generated by 
the equivalence relation relating opposite sides of the same rhombi, a 
train-track, medial cycle or thread as it can be seen graphically as a 
train-track succession of parallel edges of rhombi, graph theoretically 
as a cycle in the medial graph |11|] or knot theoretically as a thread 
in the projection of a link. 

Given the train-track of an edge (x, y) G <^ in a critical map Z, it 
defines a complex number A = 2/{Z{y) — Z{x)), of norm 2/6, which 
verifies that Exp(:A:) is null on one side of the train-track, Exp(: — 
A:) on the other. If the sum of train-tracks parallel to this one is 
non null in homology on the universal covering, then Exp(: ± A:) = 
everywhere. If it were null however, Exp(:A:) would be non trivial and 
defined independently on each connected component. 

Proposition 5.3. On a critical discrete Riemann surface, a train- 
track is never null in homology. 
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Figure 7. The image of the first sextant of the trian- 
gular/hexagonal lattice under the exponential, both dis- 
crete and continuous, for growing parameters vertically 
and finer meshing horizontally. 
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Proof [57^ . Refining by splitting each quadrilateral in four doesn't 
change neither criticality nor the homology class of a thread, so we 
can suppose that if a thread is homotopic to zero then we have a finite 
domain which completely contains it. Because all the conic angles 
are multiple of 27r, this domain is a branched covering of a region of 
the euclidean complex plane. The train-track is characterized by the 
conserved angle of each parallel side of its rhombi. As all the angles 
of the rhombi are positive, the train-track follows a directed walk and 
can not backtrack to close. Hence no train-track is null in homology. 



In the compact case, the fundamental domain is glued back into 
a closed flat surface with conic singularities multiple of 27r hence the 
fundamental group acts by a discrete subgroup G of the Galilean group 
of translations and rotations. A corollary of the proposition is that for 
all the angles appearing in this fundamental domain and their image 
by G (a discrete subset of the circle), the associated Exp(:A:) with A of 
modulus 2/6 with these given arguments must be defined as zero (or 
more precisely infinity). 

We define a region R G (}2 of a. discrete Riemann surface as convex 
iff it is connected and for every pair of faces in R along the same thread, 
if the thread is open then every face in between is also in R and if it 
is closed, the same property holds for at least one of the two halves of 
the thread. 

Given a non convex set D G (}2 contained in a disc, there is a well 
defined procedure to find its minimal convex hull by adding all the faces 
that are missing. This ultimately reduces the number of edges on the 
boundary as it doesn't introduce new threads by the Jordan theorem 
so the procedure reduces the dimension of the space of holomorphic 
functions, resp. 1-forms on it which are explicitly 00^/2 + 1, resp. 
9Do/2 — 1. Notice that 00^/2 is the number of open threads segments 
contained in D. In the tensed case (see below), a basis of holomorphic 
functions is provided by the set of functions which are non null on one 
side of a given thread segment and e. 



Conjecture 5.4. On a simply connected convex, the exponentials gen- 
erate the space of holomorphic functions. 

What is clear is that on a non convex, they don't because in gen- 
eral, not all the holomorphic functions defined on a non convex can be 
continued into a holomorphic convex on the convex hull; however. 

Proposition 5.5. In the simply connected case, every holomorphic 
function on a connected non convex which can he continued into a holo- 
morphic function on the convex hull has a unique continuation; and for 



DISCRETE PERIOD MATRICES AND RELATED TOPICS 



23 



every convex, there exists a maximal convex containing it where every 
holomorphic function on the smaller convex can be uniquely analytically 
continued. 

To prove this we need the notion of electrical moves. We could also 
prove Prop. |5]^ by implementing homotopy through Reidemeister-like 
moves with spectral parameters. They are called inversion relation 
and star-triangle relation (linked to the Yang-Baxter equation) in the 



context of integrable statistical models [13] or electrical moves in the 
context of electrical networks 



12 



6. Electrical moves 



■.P2 



II 



Pi +. P2 



III 



Figure 8. The electrical moves. 

6.1. Definition. The three types of electrical moves are presented in 
Figure ||. They can be understood as moving, splitting and merging 
the conic singularities around, keeping constant the overall curvature, 
sum on the vertices of 27r minus the conic angle. 

The first move involves the elimination of a quadrilateral with any 
discrete conformal parameter p such that two of its adjacent edges 
are identified. Let's say p labels the edge closed in a loop, pictured 
as a dashed arc. Before the move there are three vertices, of angle 
9 at the black vertex, 9' at the white vertex and 2 arctan 1/p at the 
cone summit. After the move, where the cone disappears, the angles 
are 9 — 2 arctan 1/p and 9' — 4 arctan p respectively. As 4 arctan p -|- 
4 arctan 1/p = 27r, the total curvature has not changed. In statistical 
mechanics, this move amounts to normalizing the energy by discarding 
autonomous self interaction. 

The second move replaces two quadrilaterals, with conformal pa- 
rameters pi and p2 along their respective parallel dashed diagonals, 
identified along two adjacent edges, by a single quadrilateral of confor- 
mal parameter the sum pi + p2- Dually, one replaces the series pi and 
P2 P^+p* • These formulae ring a bell to anyone versed in electrical 
networks as the parallel and series replacement of two conductances by 
a single one. The conic angles change by 2 arctan(pi + P2) —2 arctan pi — 
2 arctan P2) for the white vertices and 2 arctan(l/pj) — 2 arctan p^^p^ for 
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each black vertex. It is called inversion relation in statistical mechan- 
ics. 

The third move is the star-triangle transformation. To three quadri- 
laterals arranged in a hexagon, whose diagonals form a triangle of con- 
formal parameters pi, p2 and ps, one associates a configuration of three 
other quadrilaterals whose diagonals form a three branched star with 
conformal parameters p- verifying 

(6.1) pipi = P1P2 + P2p3 + P3pl - 



P'l + P'2+ P's 

It is the star-triangle relation, leading to the Yang-Baxter equation in 
statistical mechanics. Notice that if the central vertex is flat, then 
the equation Eq. ( |6.1|) is equal to 1 and the conic angles at either of 
the six external vertices don't change, hence flatness and criticality are 
preserved by this move. The relations for non fiat angles are more 
complicated and left to the reader. 

Given two discrete Riemann surfaces related by an electrical move, 
a holomorphic function or form defined on all the shared simplices 
uniquely defines it on the simplices which are different: 

The value of a holomorphic function 

• at the conic singularity appearing in the type I move has to be 
equal to the value at the opposite vertex in the quadrilateral; 

• at the middle point in the type II move is the < —^r—, — r— r 

^ ^ { P1+P2 ' P1+P2 J 

weighted average of the left and right values; 
at the central point in the three pointed star configuration is the 



Pi 



weighted average of the three boundary values. 



P1+P2+P3 

Moreover, the norm of the two related 1-forms (or the coboundary 
of the related functions) are equal. Hence the spaces of holomorphic 
forms for two electrically equivalent discrete Riemann surfaces are iso- 
morphic, and isometric on 1-forms. 

The notion of convexity is stable by electrical move as well. 

A discrete Riemann surface where no moves of type / or // are 
possible is called tensed [0. On a disc, electrically equivalent tensed 
discrete Riemann surfaces can be obtained from one another through a 
series of type III moves only but it is no longer the case for non trivial 
topology. Using the Dirichlet theorem and the maximum principle, one 
shows ||12[ that tensed discrete Riemann surfaces have the property 
that for every edge of the double graph A, there exists a holomorphic 
form which is not zero on that edge. It also provides with a basis of 
holomorphic functions in the case of a topological disc, together with e, 
for each open thread, there is a holomorphic function which is null on 
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one side and nowhere null on the other, unique up to a multiplicative 



constant 12 



6.2. Demonstration of Prop. |5.5| . Proof [57^. As electrical moves 
respect convexity and yield isomorphic holomorphic spaces, we can 
assume that the disc is tensed. If there were several analytic contin- 
uations then it would contradict the fact that a basis of holomorphic 
functions is given by the set of functions which are non null on one side 
of a given thread segment. The information contained in the value of a 
holomorphic function on a given convex allows to compute the value of 
the function on a greater convex. Wherever the value of the function 
is known at three vertices of a rhombi, it can be determined on the 
fourth, recursively closing the inwards corners. The result is a convex 
set where each vertex on the boundary is adjacent to at least two faces 
which are not in the set. On the square lattice = 1? for example, 
knowing the value of a holomorphic function on the horizontal and ver- 
tical axis form the origin up to (x, 0) and (0, y) allows to compute its 
values on the whole rectangle ((0, 0), (x, 0), (x, y), (0, y)). <C> 

The procedure of recursively closing the corners (and more elaborate 
conditions as well, implying more than three points) allows to try and 
extend a given holomorphic function defined on a non convex set but for 
every thread which is locally convexified, a condition for the values of 
the function has to be satisfied, if it is not satisfied then the function can 
be analytically continued into a meromorphic functions in an essentially 
non unique way. 

6.3. Continuous moves. These moves, while presented as discrete, 
are in fact continuous. A quadrilateral with a small or dually large 
conformal parameter p is associated with a thin rhombus where two 
vertices are close to one another hence the limit p = has to be un- 
derstood as an absent quadrilateral, where the edges on both sides are 
identified by pairs and dually, an infinite conformal parameter along 
an edge means that its two ends are identified into a single vertex. 

Allowing for zero and infinite conformal parameters not only enables 
to understand electrical moves as continuous, it also enables to recover 
the first and second moves as special cases of the third. Notice that 
these moves do not preserve in general the complex structure of the 
surface coded by the fiat Riemannian metric associated with the conic 
singularities. However we will see that criticality ensures that it is the 
case, the Riemann structure is unchanged by critical electrical moves. 
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Pd Pd Pl - Pl / 1 + pupp pp 

— - — - ~; — kpl-Pl) \ — 

PU Pu PR-pR PU PU 

(6.2) = i/((,'„_,,)i±4W> + 4L). 

\ Pd Pd / 

Figure 9. The continuous electrical move. 

6.4. Reidemeister moves. To recover isotopy and the usual Reide- 
meister moves, one has to allow yet another type of discrete confor- 
mal structure, namely negative parameters. Geometrically speaking, 
a negative parameter corresponds to a negative angle, or a change in 
the orientation, so the associated quadrilateral takes part in a fold, a 
overhanging wrinkle in the fabric of the surface. Notice that the change 
of all the signs corresponds to the complex conjugation at the level of 
functions or forms: a holomorphic function on a graph with negative 
parameters is a anti-holomorphic function on the same graph with pos- 
itive parameters. We will call a graph with parameters of either signs 
a virtual discrete Riemann surface. A holomorphic function on it can 
then be understood as locally holomorphic, locally anti-holomorphic 
on a non virtual discrete Riemann surface but it is not harmonic for 
this positive discrete conformal structure so we won't make use of this 
remark. Notice that the total curvature, sum on all the vertices of 27r 
minus the conic angle, is still unchanged by virtual electrical moves. 

Let's restrict ourselves to the case when all parameters on the graph F 
are equal or opposite, ±p. We associate to the planar graph with signed 
edges (F, p), a link in the tubular neighborhood of the surface IQ: The 



sign of each edge allows the resolutions of the crossings occurring in 
the thread associated to a given train-track into positive or negative 
regular projection crossings. There is of course an overall ambiguity of 
sign. The Reidemeister moves for this link are then recovered when we 
restrict the second and third moves in a certain way: 

We allow the second move only as the disappearance and appear- 
ance of a pair of quadrilaterals of opposed conformal parameters p, — p, 
leaving a collapsed quadrilateral of parameter zero. 

In terms of catastrophe theory, the whole quadrilateral with negative 
parameter is the unstable sheet of a cusp and its edges are the bifurca- 
tion lines. In this picture, the orientation of the surface is only locally 
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defined and changes across bifurcation lines wliere it is not defined at 
all. 

We restrict the third move to configurations where two parameters 
are equal and opposite to the third. 




Figure 10. The Reidemeister moves. 

Notice that the usual non virtual discrete Riemann surfaces corre- 
spond to alternating links. Alternating links are special links so not all 
the virtual discrete Riemann surfaces are Reidemeister equivalent to a 
non virtual discrete Riemann surface. 

Isotopy on the other hand is more free than the Reidemeister moves, 
the signs and spectral parameters should be irrelevant. But in fact, 
on a critical map, we can keep track of them in a consistent way. 
We no longer restrict ourselves to maps where the parameters are all 
equal or opposite and we free the restriction on the third move, but we 
keep the restriction on the second move and forbid the first move. We 
call this set of allowed moves and the associated equivalence relation 
critical isotopy. With these moves, isotopy can be performed keeping 
criticality along the way. These moves also have the property to be 
stable by refining, two critical virtual discrete conformal maps related 
by critical isotopy remain related if all their quadrilaterals are cut into 
four (or more) smaller ones. Hence, refining to the continuous limit, 
the conformal structure of a class of virtual discrete conformal maps 
(containing at least one genuine discrete conformal map) is well defined. 
Moreover, a given thread on a map yields two homotopic threads on 
the refined map obtained by splitting a quadrilateral in four. 

Notice that the first move can not occur on a critical map as it 
involves a conic singularity of a given conic angle between — vr and vr. 
Notice as well that the sum of angles at a degree three vertex, like the 
one appearing in a move of type III, can only be or ±27r on a virtual 
critical surface because it is the sum of three angles between — vr and 
TT and has to be null mod 27r by criticality. 

7. Series 



As a first attempt to define divisors, and relate them to the degrees 
of zeros and poles, we discuss here, after [O, discrete polynomials. 
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k 


3 


4 


5 


6 


7 








5 _|_ 5x'^ I 

n2 2n4 


^6 1 lOx'* 1 23x2 






1 n2 1 1 4n6 



Table 1. The first powers on the interval [0, for 
X = i/n. 



7.1. Polynomials. In a critical map Z, with Z{0) = 0, we define 
recursively the powers 

(7.1) Z'' := [ kZ^~HZ 

Jo 

with Z° = 1 and Z^ = Z. 

Consider for example <^ containing the chain {0, ,1}, the 

first Z^{x) with x = ^ are listed in Table |l| and you show in general, 
for a critical map of characteristic length 5, 

(7.2) \Z\x)-x''\<\k\x\''-^6^ 

with Afc independent of x and 5. But this constant is growing very 

rapidly with A;, namely = -j f ^if;^] where rj is the smallest 

rhombus angle in the discrete Riemann surface. It is for example not 
true that for a point close enough to the origin Z^ will tend to zero 
with growing on the contrary, if x is a neighbor of the origin with 
(0,x) G <0i and A; > 1, then Z''{x) = x'' in fact diverges with k. 
If ?/ is a next neighbor of the origin, with the rhombi (O, x, y, x') G <C>2 
having a half angle 9 at the origin, Z''{y) = ^afcU sine 00%^ e v'' 
same behavior and so has every point at a finite distance of the origin. 
It's only in the scaling limit with the proper balance given by criticality 
that one recovers the usual behavior |x| < 1 =^ {x'^l > 0. 

fc— ►oo 

See Fig. |Tl] for the comparison of the discrete and continuous powers 
on the triangular/hexagonal double. The discrepancy between the two 
increases with the degree (is null for degree up to two) and decreases 
as the mesh is refined. 

7.2. Discrete Exponential as a series. The discrete exponential 
Exp(:A:) is equal to the series whenever the latter is defined 
as its derivatives fulfills the right equation and its value at the origin 
is 1. The great difference with the continuous case is that the series 
is absolutely convergent only for bounded parameters, |A| < |. This 
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Figure 11. The image of the first sextant of the tri- 
angular/hexagonal lattice under the map z >—>■ z'' both 
discrete and continuous, for degree 1 < A; < 8 vertically 
and finer meshing horizontally. 
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suggests that asking for a product such that Exp(:A:) • Exp (://:) = 
Exp(:A + fi:) may not be the right choice. 

7.3. Ramification Number. Each rhombus F = {x,y,x',y') G <0 

is mapped by a holomorphic function / G i7°(A) to a quadrilateral 
{f (x) , f (y) , f {x') , f {y')) in the euclidean complex plane and generally 
a cycle to a polygon. If this polygon doesn't contain the origin G C, 
we define the ramification number of / at a cycle 7 to be its winding 
number around 

(7.3) 6/(7):= ' 



/(7) 



Z 



The ramification number of a face can only be { — 1,0,+!} because it 
is only four sided. If the ramification numbers at the faces inside an 
exact cycle 7 = dB are all defined, then the ramification of / at 7 is 
the sum of these numbers. The ramification number of Z at any loop 
is its winding number around O, twice for Z^. Unfortunately, things 
get more complicated for higher degrees. 

Conjecture 7.1. The ramification number of a polynomial of degree k 
is at most k, for all polynomial of degree k, there exists a cycle around 
which it has ramification number k. 

It is true at the continuous limit, it can only change by a unit, 
numerical evidence suggests that has only non negative ramification 
numbers but I don't know how to prove it in general. 

7.4. Holomorphic maps. A rhombus is mapped by a holomorphic 
function to a quadrilateral of the plane which may or may not be im- 
mersed. It is immersed whenever the images of the diagonals whether 
cross or the line going through one of them crosses the other (a kite), 
it is not immersed when the line going through them cross at a point 
outside both of the segments. Because its two diagonals are mapped 
to vectors forming an orthogonal direct basis, the only possibility for 
the ramification number to be —1 is if the quadrilateral is not im- 
mersed and the origin is inside the triangle winding the other way 



round (see Fig. 12). This case appears in the discrete analog of 1/z, in 
the derivative e', or in polynomials of degree 2 or more. 

If we don't allow this situation to occur, we can define the notion 
of holomorphic map between discrete Riemann surfaces. We need 
a more relaxed definition than criticality, namely semi-criticality 
where rhombi are replaced by (not necessarily convex) quadrilaterals. 
Note that this condition is sufficient to ensure the continuous limit 
theorem. 
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-1 



Figure 12. The three shapes of quadrilaterals and their 
ramification number. 



Let (A, p) a discrete Ricmann surface and S the associated topolog- 
ical surface. It is semi- critical for a fiat riemannian metric on S iff 
the conic singularities are among the vertices of A and every face of <) 
is realized by a linear quadrilateral flattened so that its diagonals are 
orthogonal and have a ratio of lengths governed by p. 

If the quadrilateral is convex then its diagonals are well defined as 
there exists a geodesic on the surface going from one of its four vertices 
to the opposite one, staying inside the quadrilateral, a usual linear 
segment. When the quadrilateral is in the "kite" conformation however, 
one of the diagonals would be outside the fiat quadrilateral interior and 
its geometry may be disturbed by conic singularities, it may not be a 
single segment anymore, it may even be a single point, in particular 
the notion of orthogonality may not be well defined: To compute an 
angle between two segments of geodesies, they first have to be parallel 
transported so that they cross, turning around a singularity alters this 
angle. It is why the lengths and orthogonality we are talking about 
are not those of geodesies on the surface but those computed in the 
quadrilateral isometrically flattened in the euclidean plane. 

A discrete holomorphic map f between two semi-critical discrete 
Riemann surfaces (A.p), (A',p') is a cellular complex map, orientation 
preserving, locally injective outside the vertices (where it can have a 
branching point) preserving the conformal parameters. It induces a 
genuine holomorphic map between the underlying Riemann surfaces. 
Reciprocally, a holomorphic map to a Riemann surface associated with 
a semi-critical double (A, p) from another Riemann surface such that 
the branching points are among the vertices of A induces a semi-critical 
double on the pre-image, hence a discrete holomorphic map. 

In that case, the ramification number is well defined, so that with 
B := ^/(^) the total ramification number, n = \ f~^{y)\ for any 

y in the image (counting multiplicities) the degree of the map, g and g' 
the genuses of the two Riemann surfaces, the Riemann- Hurwitz relation 
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applies and 

(7.4) g^n{g' -l) + l + B/2. 

This notion allows for the composition of discrete holomorphic maps. 

Proposition 7.2. The exponential Exp {-.X:) , for |A| < 2/6 in a critical 
map Z : U — > C, is a holomorphic map from the discrete Riemann 
surface C/ fl <02 to the semi-critical discrete Riemann surface Z{U). 

It suffices to check that every rhombus is mapped to a quadrilateral, 
whether convex or of a kite form, it will imply local injectivity as well. 
One can assume that a given rhombus under consideration contains the 
origin. By inspecting the map 

C C 

(T.5) *^ " ^ , 

one sees that the quadrilateral (1, jz^jz^, f^f ) with |2;| = < 1 
is always of that type, hence the result. 

7.5. Change of coordinate. We are now going to consider the change 
of coordinate for a critical map Z. If ( is another critical map with 
C = a (Z—b) on their common definition set, the change of map for Z^ is 
not as simple as the Leibnitz rule (a {z — b)) = Tlij=o Q) {~by 
but is a deformation of it. The problem is that pointwise product 
is not respected, Z^^^[z) ^ Z^{z) x Z^{z), in particular the first is 
holomorphic and not the second, for each partition of k into a sum 
of integers, there is a corresponding monomial of degree k. An easy 
inference shows that the result is still a polynomial in Z: 

Proposition 7.3. The powers of the translated critical map ( = a{Z — 
b) are given by 

(7.6) = a*^ ^ (''^ {-lyz'^-^B^ib) 

where (6) corresponding to V is a sum over all the degree j mono- 
mials in b, defined recursively by B^ = 1 and 

(7.7) B^{b) := (^) {-lf"''+^Z^-^{b)B^{b) 

The calculation is easier to read using Young diagrams, note the 
product of pointwise multiplication monomial 

(7.8) {Z'''{z)Y' {Z''%z)Y' . . . {Z''-{z)y" 
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B' 
B' 

B^ 



+ 1 

+ □ 

-B + 2m 

-^ + 8F + 6H-36^ + 24 mn 



+ 1-10 



+ 12 : +30 



20 ^ + 60 P + 90 ^ - 240 + 120 
'+20g- 360p + 480^ 







-90 



-90 H + 1080 - 1800 5 



rm , 



720 



Table 2. The first analogs of needed in a cliange of basis. 



with ki > k2 > ■ ■ ■ > kn as a. Young diagram Y, coding columnwise the 
partition of the integer given by the total degree k = kj x ij into 
the sum of i = Yl^=i integers. For example the following monomial 
of degree 15 = 3x2 + 2x4 + 1 is noted 

(7.9) {Z%z)y{Z^iz)yZ\z)=: ^OT. 

Then, B^{b) = c(Y)Y{b) where the sum is over all Young diagrams 
of total degree j, c(Y) is an integer coefficient that we are going to 
define and Y{b) is the pointwise product of the monomials coded by 
Y at b. The coefficient of the Young diagram Y above is given by the 
multinomials 

(7.10) c{Y) = {-l^^+' 



For example, the first Young diagrams have the coefficients ri \ i \m i i ) 



1!-1! n! 



n!, c[u 



and c(pm) 



= (-l)"+\ c(p 

("+!)! n! 



'_1 'iw+l ("+!)! 2! 
^ ^' riW. 1!1! 



(n+l)! n 



(-l)"+i2(n + 1) 
The first few terms are 



2! l!{n-l)! 2 

listed explicitly in Table 0. 

It is to be noted that the formula doesn't involve the shape of the 
graph, the integer coefficients for each partition are universal constants 
and add up to 1 in each degree. As a consequence, since Z^{z)Z^{z) = 
Z^^^{z) + 0(1/5^), 2; fixed, the usual Leibnitz rule is recovered 
in 0(1/(5^). Let's stress again that these functions B^ are discrete 
functions on the graph A which are not holomorphic. 
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The general change of basis of a given series however possible in 
theory is nevertheless complicated and the information on the conver- 
gence of the new series is difficult to obtain, even though there are some 
exceptions like the exponential, if C = (^{^ ~ (see Eq.( |5.8|) ): 

fc=0 k=0 

8. Conclusion 

The theory of discrete Riemann surfaces is shown to share a lot of 
theorems and properties with the continuous theory, extending the pre- 
vious results to the matters related to the period matrix and holonomies 
of holomorphic forms. It begs for a good definition of the order of a pole 
or a zero which would allow to prove the analog of several crucial theo- 
rems of the continuous theory, namely Riemann- Roch theorem, Abel's 
theorem and the Jacobi Inversion problem. The main challenge is to de- 
fine a good discrete analog of the exponential of a discrete holomorphic 
function so that the bilinear relations would provide Abel's correspon- 
dences between divisors and meromorphic functions. The study of the 
discrete exponential Exp(:A:Z) for a critical map Z and a constant 
A is a first step in that direction. Electrical moves are a crucial and 
powerful tool to investigate the matter. These combinatorial moves are 
clearly what is needed to make the connection with Kasteleyn theory 
of pfaffians, on the way to define discrete theta functions. 
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